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Abstract 

We observe that, at /3-deformed matrix models for the four-point conformal block, the point 
q = is the point where the three-Penner type model becomes a pair of decoupled two-Penner 
type models and where, in the planar limit, (an array of) two-cut eigenvalue distribution(s) 
coalesce into (that of) one-cut one(s). We treat the Dotsenko-Fateev multiple integral, with 
their paths under the recent discussion, as perturbed double-Selberg matrix model (at q — 0, 
it becomes a pair of Selberg integrals) to construct two kinds of generating functions for the 
g-expansion coefficients and compute some. A formula associated with the Jack polynomial is 
noted. The second Nekrasov coefficient for SU(2) with Nf = 4 is derived. A pair of Young 
diagrams appears naturally. The finite N loop equation at q = as well as its planar limit 
is solved exactly, providing a useful tool to evaluate the coefficients as those of the resolvents. 
The planar free energy in the g-expansion is computed to the lowest non-trivial order. A free 
field representation of the Nekrasov function is given. 
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1 Introduction 



Matrices have played important roles in the modern developments of quantum field theory and 
string theory. In the last several months, the /3-deformed quiver matrix models have started 
serving as a bridge [HEHSHHEIEIITIIEIEI [TUl HH H21 H31 [H] that connects the conformal blocks 
of two-dimensional conformal field theory [15] and the low energy effective actions (LEEA) of 
four-dimensional gauge theories having vanishing /3 functions. In fact, under the conjecture [T6] . 
a given Liouville conformal block gets identified with the corresponding gauge theory partition 
function dominated by the instanton sum. More than several checks have been provided, 
supporting this conjecture and its several extension [T7[IIgllIfll^l2Ill^l^l^l2^l2fll27[ 
[2H1 12HI EQl [3U [321 [331 EH EHl EH1 [371 EH]. 

These /3-deformed matrix models (or eigenvalue ensembles to call more properly), have 
Virasoro constraints at finite N as Schwinger-Dyson equations [39]. It is, therefore, not a 
surprise that the principal part of the conformal block is represented by a suitably chosen 
correlator of the model or alternatively, a potential, which has been identified with the Penner- 
type one [U SO] . 

Turning our attention to the gauge theory side, there are two faces of gauge theory LEEA 
which matrices get connected to. The one is the effective prepotential captured by the Seiberg- 
Witten curve [JT] ■ It has been demonstrated that the spectral curve of the one- matrix model 
PP, and more generally that of the A n quiver matrix model [2] are in fact isomorphic to the 
corresponding Seiberg-Witten curve written in the Witten-Gaiotto form (42j[43], giving, albeit 
a posteriori, justification of the model. The other is so-called Nekrasov partition function 
mentioned above which arises from the instanton moduli space and, under the conjecture, is 
identified with the principal part of the conformal block. The expansion coefficients of the 
former, the Nekrasov partition function [HI [45J [46] , m Qi the variable that counts the instanton 
number, ought to be identified with those of the latter in q, the cross ratio often denoted by x 
in the literature. 

In this paper, we develop computation of these g-expansion coefficients further by three 
methods that we present. The points which we have observed and will elaborate more in what 
follows are rather simple. At /3-deformed matrix models for the four-point conformal block, the 
point q = is the point where the three-Penner type model becomes a pair of decoupled two- 
Penner type models. Information on the g-expansion coefficients is well extracted at this point 
by matrix model technology and Selberg type integrals. In particular, we treat the Dotsenko- 
Fateev multiple integral [UJ as perturbed double- Selberg matrix model (at q = 0, it becomes a 
pair of Selberg integrals). We construct two kinds of generating functions for the g-expansion 
coefficients, the one for the conformal block (its principal part) and the other for the Nekrasov 
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function for £77(2) with Nf = 4, and compute some. The second Nekrasov coefficient is derived. 
A pair of Young diagrams appears naturally. 

In carrying out the computation, the integration domain of the partition function has to be 
introduced properly. In fact, the integration domain of the partition function of the conven- 
tional Penner type hermitian matrix models is over the entire eigenvalue coordinates and their 
monodromy properties are too simple to deserve a full-fledged four point conformal block la- 
belled by five generic anomalous dimensions and must be polished. The recent work of Mironov, 
Morozov and Shakirov [131 E] calls attention to the old Dotsenko-Fateev [47] multiple integral 
and the choice of paths of the screening operators. 

In the next section, we consider the Dotsenko-Fateev multiple integral and treat it as per- 
turbed double- Selberg matrix model for the expansion in q. We note an integral associated 
with the Selberg integral and the Jack polynomial to push our first method of computation. 
In sections, 3 and 4, we construct the two kinds of generating functions for the g-expansion 
mentioned above and compute some coefficients. In particular, the second Nekrasov coefficient 
is derived. The 0d-4d dictionary of the parameters is given. In section 5, computations are 
carried out from a pair of loop equations at finite iV (our second method) which are decoupled 
to each other at q = and we solve them exactly. We derive some of the g-expansion coefficients 
again as well as a few new results. In section 6, our calculation carried out at the planar level 
and the g s corrections of the loop equation (our third method) is presented. Some coefficients 
of the resolvents and the planar free energy in g-expansion are computed and are expressed in 
terms of the 4d parameters. In section 7, a free field representation of the Nekrasov function is 
given. In Appendix A, we give a summary of the 0d-4d dictionary, and some details of section 
4 are given. In Appendix B, we present computation of free energy (g-independent part) in the 
g s expansion from the Selberg integral. 

While in this paper we work on the Dotsenko-Fateev multiple integral, which is a version 
of the /9-deformed one-matrix model for SU(2) with Nf = 4, the point we have observed, 
namely, the decoupling at q = 0, holds in more general /3-deformed quiver matrix models 
[21 SU S9l [501 EU [521 [53J [5H [551 EH1 ETl [58]. This is seen particularly easily in the planar limit. 
In the planar limit, the SU(n) quiver matrix model develops a two-cuts eigenvalue distribution 
in the n = 2 case and in general an array of n — 1 two-cuts eigenvalue distributions each 
controlled by a multi-Penner potential of the following type: 

ci, a In z + c 2 , a ln(z - 1) + c 3ia ln(z - g), a = l,---,n-l. (1.1) 

The ladder structure can be seen by the form of the quiver matrix model differential Ui{z) [2]. 
While there are originally two extrema for each of the potentials, every potential at q = 
contains only a single extremum, indicating that two end-points of two cuts for each a coalesce 
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to form a one-cut distribution. See the figure for the n = 3 case. The expansion coefficients in 
q need only to be computed at an array of one-cut distributions. 



q^O, g = 2 
3 Penner 



q = 0, g = 



2 Penner 



o JJ. 



f 



f 



2 Theory of Perturbed Selberg Integral as Matrix Model 

2.1 Penner matrix model and Dotsenko-Fateev multiple integral 

The multi-Penner type "matrix model" with logarithmic potential can be defined by using a 
2d free chiral boson (f)(z) of 2d CFT with c = 1 — 6Q%, Qe = &e — (1/&e)q The multiple 
integration originates from the contour integration of N screening charges. Problem of Penner- 
type matrix model is a choice of integration contours. They may be closed or may be open. In 
[14|, taking two kinds of open contours and introducing the "filling fraction" n are suggested. 
Following the essence of their suggestion, we introduce two kinds of integration paths to define 
the three-Penner matrix model 

^3-Pcnner(g, b E ] N, Tl, «1, Oi 2 , a 3) 

:= lim g( 1 /2)(«4-2Qi5) 2 /o| . e (l/2)oi*(0) .. e (l/2)a 2 <K?) ;; e (l/2)a 3 0(l) .. e (l/2)( a4 -2Q B )<A(<?4) . 

94->-oo (2-1) 

dz : e bE<t,{z) : / dz' : e W) : |0>. 



1 We have changed our notation a little bit from our previous work [2]: &e = —ib. Normalization of the chiral 
boson has been also changed : iV2<f>(z) -> <f>(z). {(j>{z)(j>(w)) = 21og(z - w). T = (1/4) : dcf) 2 : +(Q E /2)d 2 (f>. 
With this normalization, the vertex operator : e^ 1 / 2 ) 01 ^^) ; has the scaling dimension A a = (l/4)a(a — 2Qe)- 

2 For simplicity, we consider the CFT on a sphere. Radial ordering is implicitly assumed in (|2.1[) . 
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The integer n runs from to N and serves as the "filling fraction" . Our choice of the integration 
contours C\ and C% treats the four external points 0, q, 1 and oo with high symmetry: C\ = [0, q] 
and C*2 = [1,00)3 The constant 04 is determined by the momentum conservation condition 
which comes from the zero- mode part: 

ai + a 2 + a 3 + a 4 + 2Nb E = 2Q E . (2.2) 

We will investigate the g-expansion of this four-point correlation function of CFT as a f3- 
deformed matrix model. Hence, for simplicity we assume |g| < 1. 

Let us consider the "intermediate" channel of ( 12. II) at a point z where \q\ < \z\ < 1. 
The momentum conservation at this point implies that the internal momentum is uniquely 
determined as 

«r = «i + «2 + 2nbE = — «3 — «4 — 2(N — n)bE + 2Qe- (2.3) 

Hence the only allowed intermediate states are : e^ ai / 2 ^^ ; and its descendants. By evaluating 
the expectation values of the operators, the partition function (12.11) of the /9-deformed matrix 
model turns into the Dotsenko-Fateev multiple integral 



(pq \ / ^ roc 

n / d ^ n / dz 

1=1 J o J \ J=n+1 Ji 

n 



(2.4) 



''f'-/ 1 _ - . v'f- I I ( Zj _ zi)' 2b2 '- 

1=1 l<KJ<n 
N 

x II 4 Eai {zj-q) bEa2 {zj-l) bEa3 J! {zj-zi) 2bE 

J=n+l n+l<KJ<N 
n N 

x ii n ^-^ 2b% - 

1=1 J=n+1 

Therefore, by construction, we expect that this /5-deformed matrix model (12. Ah is also a free 
field representation of the conformal block 



J r (g|c;A 1 ,A 2 ,A 3 ,A 4 ,A / ), (2.5) 



with c = 1 — QQ E an 



A i = -anion -2Q E ), (^ = 1,2,3,4), Aj = ^« 7 (a 7 - 2Q E ). (2.6) 



3 We tried several choices and judged that this choice works best. 

4 Thc vertex operator : e (i/2)(a4-2Q E )0( ?4 ) . h as t ^ e sca \[ ng dimension (l/4)(a 4 - 2Q B )(a 4 - 4Q E ), which is 
not equal to A 4 = (l/4)a 4 (a 4 — 2Qe)- Even taking into account the overall factor g^ 2 ^" 4 2 ® B > , (|2.1[) does 
not have expected scaling behaviour at finite (74. But we send 54 to 00 and in this multi-integral expression 
(|2.4[) . the parameter ol\ does not appear explicitly. The scaling dimension of (|2.4j) at 00 may be "dressed" by 
the background charge Qe- 
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2.2 perturbed double-Selberg matrix model 

In order to investigate the g-expansion of the multiple integral (12 .4p . it is convenient to make 
a change of variables 

zi = qxi, (I = 1,2, . . . ,n), z n+J = — , (J = 1, 2, . . . , N - n). (2.7) 

Vj 

The partition function Z 3 _p enner f)2.4p of the three-Penner matrix model turns into the following 
form: 



Zpert-(Selberg) 2 — <f ~ Q) 



(1/2)0203 



N L ,1 \ N L 

Yl dx I \Y[x b I Bai (l-xi) b * a *(l-qx I ) b ' ga3 ] [ \xj-xj\ 

1=1 ^° / 1=1 1<KJ<N L 

N R „i \ N R 

11 / ^ n ^ bq4 ( i - ^) b£Q3 ( i - ^^) fcBa2 n iw - vj 

7_1 </0 / r_i !/r^ 7/AF_ 



26| 



V J=1> /U / J=i 1<KJ<N R 

N L N R 

lH[(l-qx iy j) 

1=1 7=1 



Here we have renamed Z 3 _p enner , Z pert _( Selberg )2 an d introduced 

N L :=n, N R :=N-n, (2.9) 

a := -a x a 2 + 7V L + N L b E [a x + a 2 ) + Nl(Nl - 1)6|. (2.10) 

Note that, in this expression, the parameter 0:4 has reappeared through the momentum con- 
servation condition (12. 2ft : 

ai + a 2 + a 3 + a 4 + 2(N L + N R )b E = 2Q E . (2.11) 

The constrained set of the parameters (N L ,ai,a 2 ) and (Nr, a 4 , a 3 ) enters (12. 8p in a very- 
symmetric way. Indeed, if we ignore the q a part, this multiple integral (12.8p is invariant under 
the following simultaneous exchange of the parameters: 

N L i — > N R , ai < — > a 4 , a 2 < — > a 3 . (2.12) 

Note that under (T2~T2]) 

o <r-+ a' := 1 -ol a ol z + N R + N R b E {a 4 + a 3 ) + N R {N R - l)b%. (2.13) 

It is very natural to expect that the related scaling dimensions to this expected conformal 
block T{q I c ; A u A 4) Aj) are A; = (l/4)a i (a i -2Q E ) (i = 1,2,3) and A 4 = (l/4)a 4 (a 4 -2Q E ), 
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not (1/4) (04 — 2Q E )(a4 — 4Qe)- In addition to these, the scaling dimension of the intermediate 
state is given by 

A I = ^a I (a I -2Q E ), (2.14) 

where (recall fl2.3j) ) 

ai = a x + a 2 + 2N L b E = -a 4 - a 4 - 2N R b E + 2Q E . (2.15) 

Note that under the symmetry (12.121) , the internal momentum behaves as aj < — > 2Q E — aj so 
A/ is invariant under (12.121) . 

We can see that the constant cr ( I2.10p and its "dual" o' (I2.13P are equal to 

0- = Aj-Ai-A 2 , a / = A / -A 4 -A 3 . (2.16) 

Hence, the overall q factor of the partition function (I2.8P is consistently identified with that of 
the conformal block: 

q a = gAj— Ai— A2 ( 217 ) 

If we forget the Veneziano factor [59| 160] q a (l — g)^/ 2 )^^ we see ^hat a ^ g = q this integral 
decouples into two independent Selberg integrals [61]. In order to develop its g-expansion, it 
is more convenient to interpret this multiple integral as perturbation of two Selberg integrals. 
This explains why we have renamed Z 3 _ Penner , Z pert _( Selberg )2 in (12.81) . The partition function 
(12.81) of the perturbed double-Selberg model depends on the cross ratio q and seven parameters: 

Zpcrt-(Sclberg)2(g | b E ] N L , «!, OL 2 \ N R , a 4 , a 3 ) , (2.18) 

with one constraint (12.111) . 

We denote the Selberg integral by 



(N a \ n 

s N {fr,h,i)= n/ d*, n^a-^ 2 - 1 n 

\I=l J° J 1=1 1<I<J<N 

If N is a positive integer and the complex parameters above obey 



x/-xj| 27 . (2.19) 



Ref3 1 >0, Re/3 2 >0, Re 7 > -min <J i-, } , (2.20) 



N N — 1 N — 1 
then the multiple integral ( 12.19!) i s convergent [61] and equals to 

o (R o , tt r(i + irTO + (j - lhWg + (j - 1)7) , 991 , 
A ' 7) = }i r(i + 7) r(/ 3l + / 32 + (iv + J -2) 7 ) ' (2 ' 21) 
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Hence the perturbed double- Selberg model (12. 8p has a well-defined g-expansion if 

Re(b E a i )>-l, (i = 1,2,3,4), (2.22) 

/l2 s . f 1 1 Re(6 s ai) + 1 Re(6 B a 2 ) + 1 Re(b E a 3 ) + 1 Re(6 B a 4 ) + 1 

Re(6 E ) > -mm |_, -, — , — , — , — 

and \q\ < 1. 

Let us denote unperturbed "Selberg matrix model" by 



(2.23) 



-Z(Selberg)2(&E; N L , Qfj, Ct 2 ; A^r, a 4 , Qt 3 ) 

= ^selberg ip E ; N L , a u a 2 ) Z SclhcTg (b E ; N R , a 4 , a 3 ) (2.24) 
:= S Nl (1 + 6 s ai, 1 + b E a 2 ,b E ) S Nr (1 + b E a 4 , 1 + b E a 3 ,b 2 E ). 

Averaging with respect to Z (Sclbcrg )2, Zsdberg^L) and Z Se x heig (N R ) will be denoted by 
((• ■ ■ ))n l ,l r , ((■ ■ • ))n l and ((■ ■ ■ )) Nr respectively. 

Then we have an expression of the perturbed double-Selberg model: 

Zpcrt- (Selberg) 2 (<? | b E ] N L , «i, Ot 2 ] N R , Ot A) tt 3 ) 

= <f(l - g) (1/2)Q2a3 Z {Sc i bcrg ) 2 (6 £ ; iV L , a x , a 2 ; N R , a 4 , a 3 ) 

x // J](l - gx,)*— IJ(1 - qyjf^ J] IJ(1 - g*^) 26 - 

\\i=l J=l /=1J=1 II N L ,N R 

It can be rewritten as 

-^pert- (Selberg) 2 (? | A" L , «1, OC 2 \ N R , « 4 , « 3 ) 

= g A /- A i- A 2 B Q (b E ; N L , oti, a 2 ; A^, a 4 , a 3 ) #(g | AT L , oti, a 2 ; A^, a 4 , a 3 ) 



(2.25) 



(2.26) 



where 

# (&e; a i' "2; AT fl , a 4 , a 3 ) := 5^ (1 + b E ax, l + b E a 2 , b%) S NR (l + b E a A , 1 + b E a 3 , b E ), (2.27) 

B(q | 6^; iV L , ai, a 2 ; jV fl , a 4 , a 3 ) 

// Nl n r Nl n r \\ (2 28) 

= (1 _ g) (V2)a 2 a 3 // JJ (1 _ gX/)feE « 3 JJ (1 _ g!/j )to JJ _ qxiyj f<>% \\ . 

\U=1 J=l J=l J=l // Nl,Nr 

The function £?(g) = B(g | b E ; Nl, a%, a 2 \ N R , a 4 , a 3 ) has the following form of the (/-expansion: 

oo 

B(q) = l + ^2^B t . (2.29) 
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On the other hand, the conformal block of CFT has the form 
^(g|c;A 1 ,A 2 ,A 3 ,A 4 ,A / ) 

= g A '- A - A2 B { f FT \c ; A 1; A 2 , As, A 4 , A 7 ) B^\q \ c ; A 1; A 2 , A 3 , A 4) A,), 

where 

oo 

B^ FT \q) = l + J2l^f FT) - (2-31) 
1=1 

The problem is whether the following identity holds or not 

B (b E ; N L , ai , a 2 , N R , a 4 , a 3 ) = B { CFT \c ; A 1; A 2 , A 3 , A 4 , A,), (2.32) 

B(q | b E ; N L , ai , a 2 ; N R , a 4 , a 3 ) = # (CFT) (<? | c ; A ls A 2 , A 3 , A 4 , A,), (2.33) 

if we convert seven constrained parameters (b E ; N L , ai, a 2 ; N R , a^, a 3 ) into six unconstrained 
parameters (c; A 1; A 2 , A 3 , A 4 , A/) by 

c=l-6(b E -±-) , Q E = b E — (2.34) 
V oe) b E 

A l = ^a l (a i -2Q E ), (z = 1,2,3,4), Aj = i aj (aj - 2Q E ). (2.35) 

The internal momentum ai is given by (12. 151) . In the next section, we compare Be with B^ 
for £=1. 

2.3 method of calculation: from integrations to combinatorics 

We will be interested in the g-expansion of the perturbed double-Selberg model. It is a special 
case of more general perturbed Selberg model. To our surprise, the perturbed Selberg model is 
an exactly calculable model. Let us explain this for more general perturbation. Consider the 
following correlation function 

Z P ert-Seiberg(/3i, 02, T, M) ■ = S N ((3x, 2 , 7) ||exp W{xi\ #) j ^ , (2.36) 

where the averaging is with respect to the Selberg integral (I2.19p and 

00 

W(x;g) = ^2g i x i . (2.37) 

i=0 



8 



Suppose that we know the expansion of the exponential of the potential into the Jack polyno- 
mials jS21 EH rf: 

exp [^W(xr,{9i})\ = "£C?\g)P? h) (x). (2.38) 

Here P^ 1 \x) is a polynomial of x = (xi, • • ■ ,x^) and A = (Ai, A 2 , • ■ ■ ) is a partition: Ai > 
A 2 > • • • > 0. Jack polynomials are the eigenstates of 

V (xj 9 V I 7 V ( Xl + X A ( Xl 9 Xj 9 \ (2 39) 

with homogeneous degree |A| = Ai + A 2 + • • • and are normalized such that for dominance 
ordering 

P^ h \x) = m x {x) + ^a AAt m M (x). (2.40) 

Here m\(x) is the monomial symmetric polynomial. 

Let A' be the conjugate partition of A, i.e., whose diagram of partition is the transpose of 
that of A along the main diagonal. Then Macdonald-Kadell integral [671 EH [69j implies that 



,<V7>, A\ _-n(h + ( N -')i)>, ((JV+i-i)7) ; 



i>l 

X 



(A+ft + ^iV-l-ifr). 



(2.41) 



where (a) n is the Pochhammer symbol: 

(a) n = a(a+ 1) • • ■ (a + n- 1), (a) = 1. (2.42) 

Therefore, the problem of calculating the correlator (I2.36P boils down to the problem of deter- 
mining the expansion coefficients \g). 



2.4 moments in perturbed double-Selberg model 

Now, we go back to the perturbed double-Selberg model (12.251) . In this case, the parameters 
in the previous subsection are specified to 7 = b\ and 

N^N L , p 1 -H + b E a u p 2 ->l + b E a 2 , (2.43) 

5 A close connection between the Jack polynomial and the Virasoro algebra is well-known [65]. See also 
[551 155] and references therein. 
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for the "left" part. We obtain 



,«".), A\ TT ( 1 + bEai + ^ (iVl " '> 'a, C'e( N l + 1 - >) ) A , 



L=n 



( 2 + 6^(ai + a 2 ) + b 2 E {2N L -l-i)) 
- * (2 44) 

>< n 



(j . )eA (A ( -j+6|(A;-- i + l))- 

The expressions for the "right" part are obtained by replacing x = (x±, ■ ■ ■ ,xn l ) with y = 
(yi, • ■ ■ ,yN R ) and also by replacing parameters according to ( I2.12j) . 

For example, up to |A| < 2, explicit forms of the Jack symmetric polynomials P^ bls \x) 

are: 

N L 

J>2 



P$ bjs \x) = 771(1) (a;) = J2 Xl ' 

i=i 



plf E \x) = m{2) (x) + ^ mil2) (x) = ^ + ^ E (2 - 45) 

£ 7=1 s 1<KJ<N L 

PqI* s \x) = m(ia) (a;) = E 

1<I<J<N L 

Hence by using the formula ( 12. 44ft for these cases, we have 



jCV&D^XX b E N L (b E N L -Q E + ai, 
'n l ~ (a/ - 2Q S ) 



&sP ( ^V)>>„ = ^ JVM ;f JV \^^" i; , (2.46) 



h pWDwW - b E N L(b E N L + b E l )( ai + ^iV L - fe g + 26^) 
^) (b E + b E 1 )(a I -2b E + 3b E 1 )(a I -2Q E ) 



//,2 A\ _ b E N L {b E N L - b E ){ai + b E N E - Q E )(ai + b E N L -Q E - b E ) , AQ , 

2 V eP V) {x) //n l - (a I -2Q E )(a I -2Q E -b E ) " (2 ' 48) 

Here we have used the internal momenta 07 ( 12.151) to simplify expressions. 



3 Matrix Model vs Conformal Block 

In the light of the AGT conjecture, the extensive studies on g-expansion of the conformal blocks 
have been carried out in [7J [131 EEl l2Ql |27J [35] . Especially, in [131 EH] the analysis is performed 
based on the Dotsenko-Fateev integral. Therefore, we only consider the g-expansion by the (3- 
deformed matrix model ( 12.291) up to the first order. For higher order expansion coefficients, see 
these references. 
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Let us examine the g-expansion ( 12.29]) . Note that B(q) has the factorized form of the 
plethystic exponential: 



B{q) = (( exp 



°° Q k ( NL 1 \ f NL 1 ' 

- 2 E t K E x * + 2 aa 6b E^ + 2 as 

fc=i V 1=1 / V j=i / 



(3.1) 



N L ,N R 



In the expansion of this plethystic exponential, a pair of partitions (Yi,Y 2 ) naturally appears. 
The factorized form explains why in the £>-expansion 



fc=0 |Yi|=|y 2 [=A; 



(3.2) 



such that 7l(Yi) depends only on the first partition Y\ and depends only on the "left" data 
(b E , Nl) Qfi, 0^2); or equivalently (c, Ai, A 2 , A/). Similarly 7_r(Y2) is determined only by the 
second partition Y 2 and by the "right" data (b E , Nr, a 4 , a.3) or (c, A4, A3, A/). 



3.1 Hi 

From (13.11) . we can read off the first coefficient B\ as follows 



Bi = -2 I^E^+^I H^E^ + ^I • 



Using (E36D and f l2T35l) . we find 

b E N L (a x + a 2 + b E N L - Q E ) + (l/2)a 2 (ai + a 2 - 2Q E ) 



/A 1 ' 

\\ bE X/ + 2^ 2 



«3 



aj - 2Q £ 

A 7 + A 2 - Ax 
(ai-2Q B ) ' 

b E N R (a A + a 3 + b E N R - Q E ) + (l/2)a 3 (a 4 + a 3 - 2Q E ) 



Oil 



A 7 + A 3 - A 4 



(3.3) 



(3.4) 



(3.5) 



At the first order, the /3-deformed matrix model correctly reproduces the conformal block of 
CFT: 

Bi = (A / + A 2 -A 1 )(A J + A 3 -A 4 ) = ^(cft) (3 6) 

2 A/ ^ 
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4 Matrix Model vs Nekrasov Function 



In the previous section, we have considered the g-expansion of (I2.28|) which contains the 
Veneziano factor (1 — g)( 1 / 2 ) a2Q 3 as an overall factor. Instead, we can consider another q- 
expansion which does not include this factor. Let 

B(q) = (l-qY l/2)a2a3 A(q), (4.1) 

where 

II N L N R N L N R U 

A(q) := (( U(l ~ qxj) b ^ - qyj) b ^ \{ - qx lV j) 2 ^ \\ . (4.2) 

\\l=l J=l /=1J=1 I' N L ,N a 

It also has the g-expansion 

oo 

^) = 1 + J] 9 M,. (4.3) 

€=1 

This expansion is essentially the Nekrasov partition function. The analysis of the expansion of 
the Nekrasov function can be found in [2D, [22], [381 fTO] . 

4.1 0d-4d relation for the perturbed double-Selberg model 

Before going on to examine the expansion of (14.21) . we determine the 0d-4d dictionary for the 
/3-deformed matrix model. In order to obtain this, Eq. ( 14. ip at the first order in q gives strong 
restrictions on parameter relations. The first order relation of (14.11) is 

B 1 =A 1 - ^a 2 a 3 . (4.4) 



As an input, we use the explicit form of the Nekrasov function: 

[Nek , A Nek 

l (i),(o) +Ao),(i)> 



ylNek _ jNek , jNek I a r\ 



where 



iNek (a + mi)(a + m2)(a + m 3 )(a + m 4 ) 
l(1) ' (0) ~ 2a(2a + e)g* 

N( . k (a - mi )(a - m 2 )(a - m 3 )(a - m 4 ) 

" 2a(2a - ^ " (47) 



The numerator of the left-handed side of ( 14. 4p should be factorized into the product of the 
left-part data Aj + A 2 — Ai and the right-part data A/ + A3 — A 2 . This factorization can be 
explained only by the following identity: 

2a{2a + e)g 2 s 2a{2a-e)g 2 2g 2 



(2a + e)(2a-e)g 2 



a 2 - -e(mi + m 2 ) + mim 2 
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a 2 - ^e(m 3 + m 4 ) + m 3 m 4 



By comparing this relation with 



we find natural identification 

a 2 = — (mi + m 2 ), a 3 = — (m 3 + m 4 ), (4.10) 
9s 9s 

and proportionality relations 

(2a + e)(2a — e) oc A/, (4.11) 

a 2 - -e(mi + m 2 ) + mim 2 oc A/ + A 2 - Ai, (4.12) 

a 2 — -e(m 3 + m 4 ) + m 3 m 4 oc A/ + A 3 — A 4 . (4.13) 

Since A 2 = (l/4)a 2 (a 2 — 2Q E ) and A 3 = (l/4)a 3 (a 3 — 2Q E ), we can determine the 0d-4d 
relation for the /3-deformed matrix model as follows: 

a-m 2 a + m 3 

b E N L = , b E N R = , 

9s 9s 

1 / \ 1 / X 

«i = — (m 2 - m x + e), a 2 = -(m 2 + mi), 

9s 9s 

a 3 = — (m 3 + m 4 ), a 4 = — (m 3 - m 4 + e). (4.14) 

Also 0£ = 6i/g s and e = ei + e 2 , (1/&b) = — £2/ 9s- These relations convert the seven constrained 
parameters of the matrix model 

b E , N L , ai, a 2 , N R , a 4 , a 3 (4.15) 

into the six unconstrained parameters of the M = 2 SU(2) gauge theory with Nf = 4: 

ei a m 1 m 2 m 3 m 4 

— , — , — , — , — , — • (4.16) 

9s g s 9s 9s g s 9s 

Here a is the vacuum expectation value of the adjoint scalar, mi are mass parameters and ei is 
one of the Nekrasov's deformation parameter. With these parameters (14.161) . the momentum 
conservation condition (12. lip is automatically satisfied and gives no restriction on them. Under 
the exchange of parameters (I2.12p . the parameters of the gauge theory behave as follows: 

a < — > —a, mi < — > ITI4, m 2 i — > m 3 . (4-17) 

Other parameters g s , ei, e 2 , e are invariant. 
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Using the 0d-4d relation, we can rewrite the formula (12.441) in the parameters of the gauge 
theory. For example, (I2.46p . ( I2.47P and (I2.48P are rewritten as follows: 



(1/61), 



b,,pr /bB 



(1/61), 



(a — mi)(a — 1712) 
g s (2a - e) 

(a — mi)(a — m2)(fl — mi — 62X0 — m2 — £2) 

9s(ei ~ e 2 )(2a - e)(2a - e - e 2 ) 
(a — mi)(a — m 2 )(a — mi — ei)(a — m 2 — ei) 
^ 2 (2a-e)(2a-e- ei ) ' 



(4.18) 



In general, we find 



p (i/4). 



= n 

and likewise 

,(1/6|) 



(-a + m x + e x {i - 1) + e 2 (j - 1) )(-a + m 2 + ei(z - 1) + e 2 (j - 1) ) 



(4.19) 



( 2a - e - ei(i - 1) - e 2 (j - 1) ) ( ei(A$ - z + 1) - e 2 {\ - j)) 



P 



(y) 



Mr 



_1)|A| 



n 

(id)ex 



(a + m 3 + ei(j - 1) + e 2 (j - 1) ) (a + m 4 + ei(z - 1) + e 2 (j - jj 
(2a + e + ei(i - 1) + e 2 (j - 1) ) ( ei(AJ - % + 1) - e 2 (A, - j)) 



(4.20) 



4.2 expansion of A(q) 



Now let us consider the expansion of ( I4.2j) . Note that -4(g) (14. 2p takes the form of the plethystic 
exponential given by the sum of two factorized terms: 

00 „fc 



cxp 



\ / N R 

fc=l \ 7=1 / \ J=l 



(4.21) 



k=i \ 1=1 / V j=i 
Also in this case a pair of partitions (Fi,F 2 ) appears in the expansion but now A(q) has the 
form 

00 

fc=0 |Yi|+]r 2 |=fc 

Here the expansion coefficients depend on both partitions (Yx,Y 2 ) and depends on all data 
(a,m 1 ,m 2 ,m 3 ,m 4 ,g s ,e 1 ). 

For example, at the first order we have 



N L \\ // N R 
1=1 II N L W J=l //iVfl 



bE^Xl 

1=1 II A, 



«3 + b E ^2 yj 

J=X 



(4.23) 



Nr 
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We can reproduce the exact form of the Nekrasov function (14. 6 p and (14. 7p if we decompose 
(14.231) in a rather non-trivial way: 




Ao),(d = -^e^]) N+( 1+ ;^^)(( 6 *X>)) }• ( 4 - 25 ) 

See Appendix IA.3I for details. 

In the next order, the expansion coefficients A 2 can be decomposed into five terms 

A 2 = £ -4yi,y 2 = -4(2), (0) + -4(i 2 ),(o) + ^(i),(i) + Ao),(i 2 ) + Ao),(2), (4.26) 
in|+|y 2 |=2 

each of which can be written as the factorized form 



A Yl , Y2 = ^ M YuY2 (x) )) n ^ M Yu yM )) Nr (4.27) 

for some symmetric polynomials M Yl y 2 (x) and M Yl y 2 (y). At this moment, we do not have a 
clear systematic understanding on obtaining these polynomials. But at least for \Y\\ + \ Y 2 \ < 2, 
we obtained explicit forms for them. See Appendix IA.4I Suppose that one of the partition, 
say Y 2 , is empty (Y 2 = (0) ). The mass parameters and appear in -Ay^ 1 ^ through the 
following factor 

Y[ ( a + m 3 + ei(i - 1) + e 2 (j - 1) ) ( a + m 4 + e^i - 1) + e 2 (j - 1)). (4.28) 

By comparing this with (I4.20p . it agrees if A = Y\. Hence we conclude that M Yl ^(y) is propor- 
tional to the Jack polynomial P Y l J bE \y)- On the other hand, M Yl ^(x) is an inhomogeneous 
symmetric polynomial of x = (x\, . . . ,xn l ) whose highest degree is equal to |Yi|. At &e = 1 
the form of M Yl to)(%) drastically simplifies and the dominant term is given by the monomial 
symmetric polynomial: 

M Yu{0) (x) = c Y > m Y <{x) + ■■■ , at b E = 1. (4.29) 
Here Y{ is the conjugate partition of Y\ and cy is a constant. 



5 Some Computations from Loop Equation at Finite N 

In the previous sections, we used the formula (I2.4ip to calculate various objects. We show that 
some of results can be obtained by more standard matrix model technology. Using the loop 



15 



equations, the moments ((^ 7 x/)) n l and ((Y^j Uj)) n r can be determined exactly at finite N 
without using any approximation. 

5.1 loop equation at finite N 

Let us consider the loop equation at finite N of the perturbed Selberg model at the decoupling 
limit. For definiteness, we consider the left-part: 

(I! / dx i) II \xi-xj\^e^fb E ^W( Xl )\, (5.1) 

\I=1 J° ) 1<I<J<N L \ 1=1 ) 



■^Selberg (pE\ N L , Qi, Oi2, 

where 



W(x) := aii logs + a 2 log(l — x). (5.2) 



Note that here we do not introduce small g s parameter which is standard in the matrix model 
technology. The standard normalization is W (x) = (l/g s )W(x). But we do not need to consider 
g s -expansion or large Nl limit. In the next section, we write it in a more standard form which 
is suited for the ^-expansion. 
As is usual, by inserting 

1=1 1 1 

into the integrand, we obtain the loop equation at finite N: 

((^W*)) 2 )^ + (w'(z) + Qe^J (™n l (z)} Nl - fs L (z) = 0, (5.4) 

where 



Nr _ II N 



The loop equation at finite N (15.41) is valid for any potential W(z) provided that the integral 
is well-defined. The expectation value of w^ L (z) is the finite N resolvent: 



(5.6) 



For the case of our interest, the potential (15. 2p takes a very special form: double- log or 
two-Penner type. Because W'(z) = a\jz + a 2 /(z — 1), we can see that 

feWA^, (5-7) 
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where the two constants c and d are given respectively by 

Nl , \\ // v? 

C = —OL\ 



= a i w N L (0), d = q 2 // b E ^ - 1 - )) =a 2 WAr L (l). (5. 



1=1 // at l \\ /=i // 



Note that 



Z — Xt Z * — ' Z"~ 

7=1 J £=1 

where is the power sum symmetric function of xf. 

p (e) (x) = J2^i- (5-10) 

7=1 

By substituting the expansion (15. 9p into the loop equation at finite iV (15. 4p . we obtain 
infinitely many equations for the averaging of the symmetric functions of xj. At 0(1/ z), we 
have c + d = 0. At 0(1/ z 2 ), this constant is determined exactly as 

d = b E N L (a l + a 2 + b E N L -Q E ). (5.11) 

At 0(l/z 3 ), we obtain the exact expression for the average of p(i)(x): 

bENLjctj + bENL-QE) 

a l + a 2 + 2b E N L -2Q E 
which is consistent with (12.46)) . At 0(l/z 4 ) the loop equation (15. 4p at finite N contains two 
unknown moments ((b E p( 2 ))) Nh and ({(b E P(i)) 2 }) n l - 

(a 1 + a 2 + 2b E N L -3Q E )((b E p i2) (x))) NL + (((b E p {1) (x)) 2 )) NL 

(o. loj 

= b E N L (ati + b E N L - Q E ) - a 2 ({b E p {1) (x))) Nl . 

Hence beyond 0(l/z A ), the loop equation alone is not sufficient to get exact non-perturbative 
results. But it is sufficient for obtaining £>i and A±. 

We summarize the exact results obtained from the loop equation at finite iV ( 15. 4 p alone: 



{{b E p { i)(x))) NL = - ; - i oh _ M _ on _ , (5.12) 



i , V\ ///V s \\ b E N L (b E N L - Q E + oil) 

b E p m (x))) =\\b E > x i)) =7 : , , AT ^tt^> 

'I N L \\ I/ n (ai + a 2 + 2b E N L - 2Q E ) 

7 i \ b E N L (ai + a 2 + b E N L - Q E ) 
JN L {z) 



z(z-l) 

r, b E N L ( ai + a 2 + b E N L -Q E ) (5 ' 14) 
~w Nl (0) = (( b E g - ^ = , 

V ' 1 Vv & B iV L (ai + a 2 + b E N L - Q E ) 



^(i) = ((&*Ei3 



/=i V/jVx, 



Except the first equation, which is (12.46)) . the remainder is a set of independent non-perturbative 
results for the perturbed Selberg model. 
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6 Computation at the Planar Level and g s Corrections 

In this section, we consider the (^-expansion of the loop equation (15.41) . For this purpose, we 
introduce the scale g s and rescale various parameters. First, we set 

oti := g s aa, i = 1,2,3,4. (6.1) 

Then, the potential is rewritten as 

W(x) = —W(x), W(x) := fiilogx + 5 2 log(l -x). (6.2) 
9s 

Now the perturbed Selberg model takes the standard form in matrix model calculation: 

n f dx i) n \*i - x j\ 2bl ex p ( — e w ^n ■ (°) 

K I=l^° J 1<I<J<N L \9 S 1=1 J 

A similar expression holds for the right part. 

We will use the 't Hooft coupling Sl and Sr and their deformed version: 

S L := g s N L , S R := g s N R , S L := g s b E N L , S R := g s b E N R . (6.4) 

In the rescaled parameters, the momentum conservation condition (12. lip turns into 

(k\ + a 2 + a 3 + a 4 + 2S L + 2S R = 2e. (6.5) 

We also change the normalization of the resolvent wn l (z) := g s w NL (z), 

w Nl {z) = (^b E g s J2jz^jJ > WN *( z y- = Y E 9'^J^jj ■ ^ 6 ' 6 ) 
Also, we set 

f NL (x) :=gJ NL {x) = (^b E g s 2_^ — — Jj . (6.7) 

Now the loop equation at finite N takes a more standard form: 

)) +W'(z)w NL (z) + ew' NL (z)-f NL (z) = 0. (6.8) 

7 N L 

Here e = g s Q E — ei + £%■ A similar expression holds for the right part. 
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6.1 b E = l (e = 0) 

We first consider = 1 case. In the planar limit, ( 16. 8 j) reduces to an algebraic equation 

^lW 2 + W'(z) - /(*) = 0, (6.9) 

where 

w L (z) = lim w Nl (z), f(z) = lim f NL (z). (6.10) 

N L -+oo N L —>-oo 

In the planar limit, the momentum conservation condition (I6.5P becomes 

a 1 + a 2 + a 3 + a 4 + 2S L + 2S R = 0. (6.11) 
The solution to the planar loop equation (I6.9P is 



W'(z) IfW'iz^ 2 



w L (z) = ±J. + J\—±J-J +f( z ). (6.12) 

The form of f(z) is already restricted to 717 + 7 and the residue at infinity, which is proportional 
to c + c', vanishes due to the stationary condition on the distribution of the eigenvalues. Hence 

m = (6.13) 

Picking the residue of wl(z) at infinity, and equating it with Si = g s N L , we obtain 

c = Sl + (a 1 + a 2 )S L . (6.14) 

The solution takes the form of 



1 ( ol\ ai 2 \ ll f an a 2 \ 2 , S L {ai + a 2 + S L ) 
w L {z) = w(z; an, ac 2 , S L ) = 1 + \ h 



(6.15) 

and is completely determined. Likewise, 



- c , , 1 f 014 a 3 \ llfai 4 . «3 V , S r(^ + «3 + S R ) 
w r {z) = w(z;a4,a 3 ,b R ) = 1 +4/- 1 + 



2 \ z z — 1 / V 4 \ z z — 1 J z[z — 1 J 

(6.16) 

We expand them as follows 

00 (f) 00 {£) 

fcO £=0 
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They correspond to 

At the first order, the coefficient of (14.211) and that of (13.11) are respectively 

A 1 = -~ 2 (wf + 6 2 )w^ - - 2 wf{w^ + o 3 ), (6.19) 

Si = -2 + ^2) Ug } + i« 3 ^ = A - ^«2« 3 - (6.20) 
For simplicity, we use 

a/ := ai + a.2 + 25*/, = —64 — 63 — 25l. (6.21) 

At the first order, we obtain 

m = S L (a 1 + S L ) m = _g fl (a 4 + fl0_ (6>22) 
a/ a/ 

Thus 

(i) . ~ («2 + S L )(ai - S L ) (i) _ (a 3 + 5 fl )(5/ + 5 fl ) , . 

+a 2 = z , u>V + a 3 = z , (6.23) 

oj (Xi 

and 

£fl(d 4 + 5 , fi )(a 2 + S l )[oli - S L ) S L (&i + S L )(6 3 + S R )(aj + S R ) 

At = 2^2 • ( 6 - 24 ) 

9% oq 9s a i 

Adding — (1/2)6263/ g 2 s coming from the Veneziano factor g( 1 / 2 ) Qia2 (l — ^(v 2 )^^ we obtain 

Bi = ^-^2 (^Sl(6i + a 2 + S L ) + 62(61 + 6 2 ) ) (2S R (6 4 + 6 3 + S R ) + 63(64. + S3)) . (6.25) 
In the planar limit, the higher coefficients At and Be (I > 2) are factorized as follows: 

At = B t = ±(B0'- (6-26) 

The planar contributions to At and exponentiate to give 
00 00 

^planar (?) = ^ <f At = exp (g^) , -B p i anar (g) = ^ = exp (gBx) . (6.27) 

£=0 £=0 

The higher moments (16.181) in the planar limit are, therefore, not really relevant for the calcu- 
lation of At and Bt (I > 2), but are useful for the calculation of the planar free energy. For 
simplicity, let us consider 

A(q) = exp (T(q)) = exp (^o(<?) + Fi(q) + ^ 2 (?) + • • • V (6.28) 
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Then, we can easily evaluate the g-expansion coefficients of the planar free energy Fo{q) 

oo 

^o(g) = (6.29) 

k=l 

We already have 

2 + S L )(a 3 + S R )(a>i - S R ) {S L + a 2 ){ai - S L )S R (a A + S R ) 
•n>,i = S's^i = =5 T2 • (6.30) 



From (I4.2ip . we easily see that 



J 7 Q,k 



a 



a 



^ + a 2 )ttV + W L ( W R + a V 



In the next order 



w 



(2) 



61,(01 + S , L){ai(«i + a 2 ) + Sl(3«i + 2a 2 + 3Sl)} 



^,3 



a 



(6.31) 



(6.32) 



(2) 

and a similar expression for w R is obtained by the replacements 



S L — > S R , &i — > &4, a 2 — > &3, 



—ft/. 



Then 



(2) . ~ («2 + S L )(ai - S L )[ai{pti - S L ) + S L (a 2 + S L )} 
w K L ' + a 2 = - — 



a. 



(2) 

and a similar expression is obtained for w R + a 3 . By substituting these relations into 



0,2 



>[ 2) + « 2 )^ 2) + 4 2) (4 2) + «3) 



(6.33) 



(6.34) 



(6.35) 



we can obtain an explicit form of ^0,2- 

In this way, w^ 1 and for k = 1,2,3,... can be computed and upon substitution 
into f!6.3ip . J"o,fc are determined. The g-expansion coefficients of the planar free energy are 
completely controlled by the single algebraic function w(z; a, (3, S). 

We can match parameters and can get planar 0d-4d relations also by using the planar loop 
equations. For example, let us match parameters of matrix model with those of CFT. In the 
planar limit, we have an expression of B\ (16.251) . On the other hand, 



(cft) _ (A/ + A 2 - Ai)(A/ + A 3 - A 4 ) 



2A, 



(6.36) 



We know that in the c = 1 free boson system, the vertex operator : Q^ a i/ 2 9s)<t>( z ) . ^he scaling 
dimension Aj = 1/(4^)5?. By identifying the internal momenta with ft/ = (l/g s )6ci, f l6.25|) 
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completely agrees with the CFT expression (I6.36P @. As for the gauge theory parameters, 
analysis can be done in the same way as that of section HI The planar 0d-4d relation is 

S L = a-m 2 , S R = -(a + m 3 ), 

oil = m 2 — mi, a 2 = wi2 + mi, 

a 3 = m 3 + m 4 , a 4 = m 3 — m 4 . (6.37) 

Also, a i = 2a. Using the relation (16.371) . the solution (I6.15P to the planar loop equation takes 
the form 



1 (mi — vni m>2 + mi\ 1 (vri2 — m\ 1112 + mi\ 2 (a + 1112) (a — 1712) 

WLix) = -2{—— + ^r) + \Ji[—— + ^r) + — 4^1) — • 

(6.38) 

A similar expression for wr(x) can be obtained by the replacement (I4.17p . 
First few moments (16.221) . (16. 23 p . (I6.32f) and (16.341) are rewritten as follows 





(a — rrii)(a — 1712) 
2a 










(1) , ~ 
w y L + a 2 = 


(a + m 1 )(a + m 2 ) 
2a 










(2) 


(a — mi)(a — m 2 ) 


(a + mi){a 


+ m 2 )(a 


— mi){a 


- m 2 ) 


2a 




8a 3 






(2) , ~ 

w y L + a 2 — 


(a + mi)(a + m 2 ) 


(a + mi)(a 


+ m 2 )(a 


— mi)(a 


-m 2 ) 


2a 




8a 3 







(6.39) 



and similar expressions for the right part. 

First two of the planar instanton contribution (I6.3ip obtained from the expansion coefficients 
of fhMD are 

_ (a + m x )(a + m 2 )(a + m 3 )(a + m 4 ) (a - m x )(a - m 2 )(a - m 3 )(a - m 4 ) 

Ai = ^ + ^ > ( 6 - 40 ) 



r rii( Q + m A (a - mjja - m 2 ) \ ( 
• Fo ' 2 ~ 8^ I 1 4^ J ^~ 



(a — m 3 )(a — m 4 ) 



4(7 2 

7 (6.41) 

Yli(a-mi) I (a + m 1 )(a + m 2 )\ f (a + m 3 )(a + 



V- V ' 4a 2 y V I" 2 

6.2 6^ 7^ 1: ^-corrected resolvent 

We briefly comment on the case of 6e 7^ 1 and on the g s -correctior]ll to the resolvent ( 16.61) . We 
choose g s and as independent parameters and consider the perturbation of the loop equation 



6 If g s — > while keeping 04 finite, the external vertex operators carry very large momenta ai — (1/ ' g s )a.i. 
7 It is known that [70j when bs 7^ 1, the ^-expansion contains odd powers in g s . 
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( 16.81) in g s . The planar results for &g 7^ 1 can be obtained from those of &e = 1 by replacing 
the 't Hooft parameters Sl, Sr by the deformed 't Hooft parameters Sl, Sr respectively [2]. 
(See d5aj). 

In the g s -expansion, some care is necessary on the choice of parameters to be fixed. We 
would like to compare the g s -expansion coefficients with those of the gauge theory. We have an 
exact 0d-4d relation (14.141) and we fix the parameters of M = 2 gauge theory. There is a subtle 
point. Notice that for frg 7^ 1 (i.e., Qe 7^ 0), the parameters in the matrix model potential 
(16.21) are rewritten as 

&i = m 2 - mi + g s QE, a 2 = m 2 + mi, 

a 3 = m 3 + m 4 , a 4 = m 3 — m 4 + g s QE- (6.42) 

The parameters ot\ and <5 4 are inhomogeneous with respect to the degree in g s . The derivative 
of the potential has the following ^-expansion 

W(z) = K(z) + 9 -^, Wi(z):=^^ + ^±^.. (6.43) 

For simplicity, we consider the left part only. 

Recall that we have the exact form of the function Jn l { z ) (16.71) (see (I5.14p ). We see that it 
receives no g s -correction: 

, , . (a + m2)(a-m 2 ) . , s ln ... 

In l {z) = ~ ,J_ ^ " = h(z). (6.44) 

The resolvent (16.61) admits the following g s -expansion 

w Nl (z) = w L (z) + g.w£ /2 \z) + g 2 s w { l\z) + ■■■. (6.45) 

Up to 0(g s ), the first term in the loop equation (I6.8P behaves as 

(({w NL (z)) 2 )) NL =w L (zf + 0(g 2 s ). (6.46) 

By substituting ( I6.43p . ( I6.44|) . (I6.45P and (I6.46P into the loop equation ( I6.8p . we can easily see 
that the 0(g s ) relation is given by 

(2w L (z) + Wl(z))w% /2 \z) + ^w L (z) + Q E w' l (z) = 0. (6.47) 
Hence the 0(g s ) corrected part is easily determined as 

»^M = - «(2M»?W)) 5 ( ™ lW) - (6 ' 48) 

Therefore, the ^-correction is also characterized by the algebraic function w(z; a, (3, S). Here 
wl(z) takes exactly the same form as the planar solution (I6.38P to 6g = 1 case. 
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7 Free Field Representation of Nekrasov Function 

Finally, we come back to the definition of the perturbed double-Selberg model (12.11) given in 
terms of the free chiral boson. Using this definition, we can write a free field representation of 
the Nekrasov function. We write the modes of the free boson as 



■o + a log z - Y] — z n = 4> + a logz + <f> + (z) + <j>-(z), (7.1) 



Z — < n Z — < n 

n>0 n>0 

The non-trivial commutation relations are (cf. footnote 1) 

[«o, M = 2, [a n , a m } = 2n 5 n+mfl . (7.3) 

We decompose the Virasoro generator Lq into the non-zero mode part and the zero-mode 
part: 

L = K + A, £:=l^a_ n a n , A := ^a 2 - ^ a . (7.4) 

n>l 

The eigenstate of cto is defined by 

|a):= e ( a / 2 )*>|0), («| := (0|e- (a/2) ^. (7.5) 
Let us introduce "left" intermediate state \ai',q) and "right" intermediate state (ot/| by 

| a z ; q) =: e (1/2)Q ^ (<?) : (Q[o, q] f L | «i>, (a r \ = (2Q E - a A \ (Q {1 ,oo]) Nr : e^" 3 ^ :, (7.6) 



where 



Note that 



QWM = f" d* : J*« M) : . (7.7) 
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a \ai;q) = \ar,q) a/, a/ = a x + a 2 + 2b E N L , (7.8) 

A| af, q) = | a/; g) A/, A 7 = -a/(a/ - 2Q B ). (7.9) 

The definition of the perturbed double-Selberg/three-Penner matrix model (12. ip can be con- 
cisely written as 

Zpert-(Selberg) 2 = («/ !«/;<?). (7-10) 

The dependence on the cross ratio g appears only in the left intermediate state. 
Recall that the vertex operator transforms under the scale transformation 

q Lo (: e (a/2)0(2) :)q~ L ° = q Aa : e^ 2 ^^ : . (7.11) 
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The "left" screening charge Q[o, q ] behaves under the coordinate transformation z = qx as 

Qlo, q] = Q Lo U dx : e 6 *^ q~ L ° = q L ° Q m q~ L °. (7.12) 



Hence the left intermediate state can be rewritten as 

\ ai ) = q L ° : e VW**0-) : (Q [0jl] )^ | a x ) q^^ 

= q L °- Al ~ A2 \N L ,a 1 ,a 2 ;b E ), 

where 

I N L - ai , a 2] b E > =: e^WW . (Q [0>1] )^ | ai) 

'N L i \ jV L 

nf** n 

v^l" 70 / j=l 1<KJ<N L 

x exp 



(7.13) 



J] / dx, IJrc^a-arj)^ J] Izj-zjI 

' ' " ' ' 1 1<KJ<N L 

(l NL \ 
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(7.14) 



We call |iV L , ai, a 2 ; &e) Selberg state. Similarly we can introduce the "right" Selberg state 
(N R ,a 4 ,a 3 ;b E \. 

Then the partition function of the perturbed double-Selberg model can be compactly written 

as 

Zp Cr t-(seiberg)2(<?) = <? Al ~ Al ~ A2 (N R , a 4 , a 3 ; b E | q K | N L , a x , a 2 ; b E ) • (7.15) 

Furthermore, if we convert matrix model parameters into the those of gauge theory by the 
Od — 4d relations, we also rewrite it as follows: 



Z 



2 



pert— (Selberg) 

= q a ( | e ( 1 / 2 9s)( m 3-m.4-e)0o ^ -(l/&E3 s )(a+m 3 ) . e (l/2g s )(m 3 +m 4 )</>(l) . (7-16) 

X q^ : e ( 1 / 2 9s)(m2+m 1 )(/,{l) . \(l/6 B p,)(o-m a ) e (l/2ff s )(m 2 -mi+e)0o | Q y 

This gives a /3-deformed version of the formula seen in [59] . 
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A Matrix Model and Nekrasov Partition Function 



A.l 0d-4d Dictionary 

In this appendix, we summarize the 0d-4d dictionary for the /3-deformed matrix model (the 
perturbed double-Selberg model) (12. 8p . The following is our notation for parameters of Od 
matrix model and 4d gauge theory: 



^pert-CSelberg)^ I ^ N L, «1, Ot 2 ] N R , « 4 , « 3 ) = ^Nck (q 



£i a mi m 2 m 3 m 4 

9s 9s 9s' g s ' 9s' g s 



(A.l) 



In addition, on 4d side, the following symbols are often used: 

-- = -1, e = e 1 + e 2 . (A.2) 

9s 9 s 

The matrix model parameters obey the momentum conservation condition: 

ai + a 2 + a 3 + a 4 + 2(N L + N R )b E = 2Q E . (A.3) 

Also, on Od side, we use 

Q E = b E - —, ai = ai + a 2 + 2b E N L = -a 3 -a 4 - 2b E N R + 2Q E . (A.4) 
o E 

The 0d-4d relation is given by 

a-m 2 a + m 3 

b E N L = , b E N R = , 

9s 9s 

ai = — (m 2 - mi + e), a 2 = -(m 2 + mi), 

9s 9s 

a 3 = — (m 3 + m 4 ), a 4 = — (m 3 - m 4 + e). (A. 5) 

9s 9s 

A.2 some relations among parameters 

Some relations among parameters are collected below: 

— = b E , — = -T-, — = = Qe = b E - — , (A.6) 

9s 9s b E g s g s b E 

11 1 

- = o(<*j - Qe) = + a 2 + 2b E N L - Q E ) = ~(-a 3 - a 4 - 2b E N R + Q E ), (A.7) 
g s Z Z z 

— = -z{oi2 -cxi + Qe), — = -( a 2 + «i - Qe), 
9s 2 g> s 2 

m 3 1 m 4 1 

— = ^(«3 + «4 - Qe), — = ^(a 3 - «4 + Qe)- (A. 8) 
9s z g s Z 
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a + mi a — m 1 

= a 2 + b E N L , = ai + b E N L - Q E , 

9s 9s 

a + m 2 a — m 2 

= ai + a 2 + o B 7V L - Q S) = b E N L , 

9s g s 

a + m 3 a — m 3 

= -b E N R , = -a 3 - a 4 - OeA/r + Q E 

9s 9s 
a + m 4 



-a 4 - 6^ + Qe, = -« 3 - b E N R . (A.9) 

a/ = ai + a 2 + 2b E N L = -a 3 — a 4 — 2b E N R + 2<5s, 



(A.IO) 



2a - e 

= a/ - 2<5e = «i + a 2 + 2b E N L - 2Q E = -a 3 - a 4 - 2b E N R . 

9s 

A.3 equivalence of dOD and (H25D with (USD and ( ICT ) 

If we use the following relations, 

1 _ aj - 2Q g _ 2a - e 1 | Qe _ ai _ 2a + e (All) 
oci — Qe oii — Qe 2a a/ - Q E ai~Q E 2a 

h \ - b EN L (b E N L -Q E + ai) _ (a - mi)(a - m 2 ) ,^ , 

E h Xl // NL ~ (<*i-1Qe) ~ 9s(2a-e) ' 1 " j 

^ \\ b E N R (-a4-b E N R + Q E ) (a + m 3 )(a + m 4 ) ,„> 

°» + (^)f|'l.- ( ' +B S +w,) | (A ' 14) 

* + (^)fS^-- &z2 S z ^ (A15) 

we can see that (I4.24p and (I4.25P reproduce the Nekrasov function (14. 6 p and (14.71) respectively. 



27 



A.4 A 2 

The second order g-expansion coefficient A2 of (14. 2 j) is given by 

\ 2 \\ // / ^fl X 2X 



2 

1 

2 



2/J 



1=1 / // Ar L \\ \ J=l / //jV^ 

&b E xj ) » (( I a 3 + b E ^2 VJ 

1=1 J II N L W V J=l J II N R 



1 

~ 2 
1 

~ 2 

By using 



J=i / \ /'=i / // n l \\ V J=i /V J'=i / // JV H 

JV£ \\ // N R 

b E ^xj + a 2 J) (I b E Vj 

1=1 II N L \\ J=l //jVfl 

A^L \\ // iVfl 

/=! //jv \\ J=l //jv 



= W)(*) = ^ - rf|| p ( ( i) 6|) ^)) > ( A - 18 ) 

and similar relations for the right part, we can convert the symmetric polynomials in (1A. 16j) into 
the Jack polynomials. Through (|4.18|) . the exact expressions for these moments are obtained: 



N L t Utl - t2jl ^-t n ^- t - tlj 



(a — mi)(a — 


m 2 )(a - 


m l — ~ m 2 — 




ei(ei - 


-e 2 )(2a- 


-e)(2a-e-ei) 




(a — mi)(a — 


m 2 )(a - 


mi - e 2 )(a - m 2 - 




62(61 - 


e 2 )(2a- 


■e)(2a-e-e 2 ) 




(a — mi)(a 


- m 2 )(a 


— mi — ei)(a — m 2 




^(ei 


-e 2 )(2a 


-e)(2a-e-ei) 




^ (a — mi)(a 


- m 2 )(a 


-mi - e 2 )(a - m 2 


-62) 



7=1 ll Nl ^ ^ (A.20) 



By substituting these formula into ( 1A.16I) . the expansion coefficient *4. 2 is obtained. Since the 
result is rather lengthy, we do not write it here. But we have checked that it exactly coincides 
with the sum of the five Nekrasov functions: 

A2 = -4(2 G ) k (o) + -4.(12), (0) + -4(i) k (i) + -4(o) k (i 2 ) + -4(o C ) k (2)- (A. 21) 

Nek 



An explicit form of these Ay*y 2 can be found, for example, in [T8] . 
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A. 4.1 qualitative explanation on rearrangements 

Unfortunately, we do not have a systematic algorithm for rearrangements of the exact expression 
into the Nekrasov functions (14.261) such that each term has a factorized form (14.271) . But a 
general pattern can be easily read out from flA.161) . Note that the mass dependence of the 
moments flA~T9|) and fTA~20l) are given by certain symmetric polynomial of a — mi and a — m 2 . 
They do not contain a + mi and a + m 2 . Likewise, on the right part, the moments depend on 
a + rris and a + symmetrically. The sign flip of mass is caused by mixture with the constants 
a 2 = (rrii + m 2)/g s (or a 3 = (m 3 + ra^)jg s for the right part). For example, let us consider the 
first term in flA.16j) : 



(A.22) 




The term on the right part contains the Jack polynomial PLj E (y): 




b 2 E P$ b% \y) + ■ ■ ■ , (A.23) 



thus this "dominant" part gives 



b E x {{b E P^\y))\ =- {a + m3Ka + ^ + ™ 3 , + £2)(a + m4 + £2) . (A.24) 
W (2) yyi l'N R e 2 (ei -e 2 )(2a + e)(2a + e + e 2 ) 



In the left part, a quite non-trivial portion M( 2 ) ( )(x) of the total terms gives a "dominant" 
contribution: 

I • 1 L | / \ L 

J- 9 . , \ . 1 / , \ 

Xi 



V 1=1 J 1=1 V 1=1 



(A.25) 



^(2),(o)( a; ) + sub-dominant terms, 



where 



\* / \ 1 2 "2 . (2a-ei)a 2 / D (i/&|), x 
M ( 2),(o)(x) := -a 2 - ^- + (2a + ^ (^(i) (*, 

(2a " 6 ~ e2) ^|P^(x)) (A.26) 



e 2 (ei-e 2 )(2a + e 2 ) V E d 2 ) 
( ei + e 2 )(2a - e x + e 2 )(2a - e - eg) / p (i/b%) ( 
f 4a^ s (2a + 6 2 ^ V (2) W 



Here, we have determined this polynomial by requiring the linear factorization of the moment 
(as in 0A.24P ) and also by the assumption on mass dependence only on a + mi and a + m 2 
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symmetrically. Indeed, while the form of the polynomial flA.26j) is messy, its moment takes a 
quite simple form: 

/ x\\ (a + m 1 )(a + m 2 )(a + m 1 + e 2 )(a + m 2 + e 2 ) , . 

M (2),(.Q)\ X ) /) = a 7^ — , \ • (A.Z7J 

'I n l 4eie 2 a(2a + e 2 ) 

(l/b 2 ) 

Hence, seen at the level of moment it is essentially dual of the Jack polynomial P,^ E (x) such 
that the sign of a is reversed. 
Now by setting 

M m o)(y)--=b 2 E P^ bl \y), (A.28) 
we can see that the first term ( 1A.22I) contains the portion 

M (2)(0)( x)\\ (m {2) m)) = nU fl+ "?X fl +7 +e2 > (A29) 

( 2 ),(0)v J// Nl \ {2),(0)^J// Nr 4 eie 2 a (2a + e 2 )(2a + e)(2a + e + e 2 )' V ; 

which is exactly equal to -4.^w y 

Note that at b E = 1, the symmetric polynomial M^),{q){x) (1A.26j) turns into the form 



M mo) (x) = ^a 2 2 - ^a 2 + a 2 P${x) + P^{x). (A.30) 

It has quite simple coefficients and the dominant term is P$Jx) = m^(x), the monomial 
symmetric polynomial characterized by the conjugate partition Y[ = (l 2 ) of the partition 

Yi = (2)- 

A. 4. 2 polynomials M Yll Y 2 (x) and M Yl y 2 (y) 

Here we write the polynomials My 1 y 2 (x) and their pairs My 1 y 2 {y) (I4.27P with |Yi| + \ Y 2 \ = 2. 

For Yy = (2) and Y 2 = (0), they are given by (IA.26I) and (1A.28|) respectively. While for 
Y x = (l 2 ) and Y 2 = (0), 



M (ia) , (0) (aO = \a\ + h E a 2 + (^f) c^b E P^'\ 



e(2a + ei - e 2 )(2a - e - e x )\ &2 p (i/&|), . 



2( ei -e 2 )a(2a + ei ) ^ £ ( l2 ) 
• 2a + ei J^ 2 ) 



(A.31) 



Also, for Fx = (1) and Y 2 = (1), 

>(i/b|)^ ^e 2 (2o-e-ei) nh2 D (i/6|) 
(ei - e 2 )(2a - e x 



^ l 0(1/4)/ \ 9s^2{^a - e - ei 2 , 

^(D,(i)(z) = a 2 g s b E P {1) (x) —————— —2b E P (ia) 



2a-e-e 2 \ p (i/&|) M 
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(A.32) 



The rest can be obtained easily from these formula by the exchange My li y 2 (a;) < — > My 2 ,y 1 {u) 
along with the replacement of parameters (I4.17P . 
The moments for My 1 y 2 (x) are given by 

,, , N \\ (a + mi) (a + m 2 )(a + m 1 + ex) (a + m 2 + e x ) 
M mm ( X ))) Ni = 4 9 | o( 2a + ei ) ' 

(a + mi) (a + m 2 )(a - mi) (a - m 2 ) 



n l g s (2a - ei)(2a - e 2 ) 

Note that at = 1, 

M (l 2 ),(0)(^) = + 2 "2 + «2P ( ( 1) ) (X) + Pgl 



(A.34) 



M(i) j(1 )(z) = g s a 2 Pfi(x) + g s Pfy(x) + g s P^(x). 

B Free Energy from Selberg Integral 

In this section, we briefly discuss the contribution to the fee energy of the perturbed double- 
Selberg model f)2.25p . i.e., -Z^seiberg) 2 : 

^(Selberg) 2 = exp (T L + T R ) , (B.l) 

where 

J r L = \ogS NL (l + b E a 1 ,l + b E a 2 ,b 2 E ), T R = log S Nr (1 + b E a 4 , 1 + b E a 3 , b%). (B.2) 

Since we have an exact expression of the Selberg integral Sn(Pi, P2, 7) (12. 21 j) . an integral rep- 
resentation of the free energy is easily obtained: 



(B.4) 



where 

Q L (t) = exp (-b 2 E t) + exp (-b E acit) + exp (-b E a 2 t) 
- exp(— (1 + b E (ai + a 2 ) + (N L - l)b%}tj, 

and similar expression for Tr. 

B.l ^-expansion of the free energy 

For simplicity, we assume the deformation parameter b E real and positive. Then, the free energy 
( IB. 31) can be rewritten by using the parameters of the gauge theory. Through the relation (14. 141) . 
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we have 



COO Jy. 







t 



(1 - e-( a - m2 ) 4 )^ L (t) (a - m 2 ) ( (1 + e - * 6 **) ;| 



(B.5) 
where 

Q L (t) = e~ 9sbEt + e ~^ m2 ~ mi+ 9sQE)t _j_ e -(m 2 +mi)i _ e -(a+m 2 )t g-j 

Let us consider the g s -expansion of ( IB. 51) : 

= £ ^" 2 4 fc/2) = + 1 4 1/2) + 4 1} + • • • • (B.7) 

fe=0 ^ s & 

For example, the planar part is given by 

^i 0) = -E|(a + mi) 2 log(a + mj) + (a - m;) 2 log (a - m;)j 



8=1 



— 2a 2 log(2a) — - (m 2 + mi) 2 log(m2 + mi) — - (m 2 — mi) 2 log(m2 — mi). 
r(o) 



(B.8) 



Combining with F R , we obtain a perturbative part of the prepotential for the gauge theory. 
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